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Abstract 

According to quantum chromodynamics, matter at ultra-high density and low temperature is a 
quark liquid, with a condensate of Cooper pairs of quarks near the Fermi surface ("color super- 
conductivity")- This paper reviews the physics of color superconductivity, and discusses some of 
the proposed signatures by which we might detect quark matter in neutron stars. 



1. Introduction 

One of the most striking features of quantum chromodynamics (QCD) is asymptotic freedom: 
the force between quarks becomes arbitrarily weak as the characteristic momentum scale of their 
interaction grows larger This immediately suggests that at sufficiently high densities and low 
temperatures, matter will consist of a Fermi sea of essentially free quarks, whose behavior is 
dominated by the high-momentum quarks that live at the Fermi surface. 

However, over the last decade it has become clear that the phase diagram of dense matter 
is much richer than this. In addition to the hadronic phase with which we are familiar and the 
quark gluon plasma (QGP) that exists at temperatures above about 170 MeV, there is a whole 
family of "color superconducting" phases, which are expected to occur at high density and low 
temperature; for a detailed review, see Ref. [1]. These phases have observational importance, 
because they may occur naturally in the universe, in the cold dense cores of compact ("neutron") 
stars, where densities are above nuclear density, and temperatures are of the order of tens of keV. 
It might conceivably be possible to create them in future low-energy heavy ion colliders, such 
as the Compressed Baryonic Matter facility at GSI Darmstadt [2]. Up to now, most work on 
signatures has focussed on properties of color superconducting quark matter that would affect 
observable features of compact stars, and we will discuss some of these below. 

2. Color superconductivity 

2.1. Cooper pairing of quarks 

The essential physics of color superconductivity is the same as that underlying conventional 
superconductivity in metals, and also superfluidity in liquid Helium, nuclear matter, and cold 
atomic gases. The crucial ingredients are a Fermi surface and an attractive interaction between 
the fermions. Quark matter has exactly these ingredients. It was shown by Bardeen, Cooper, 
and Schrieff'er (BCS) [3] that if there is any channel in which the fermion-fermion interaction is 
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attractive, then there is a state of lower free energy than a simple Fermi surface. That state arises 
from a compUcated coherent superposition of pairs of particles (and holes) — "Cooper pairs". 
This can easily be understood in an intuitive way. The free energy at zero temperature is F = 
E-juN, where E is the total energy of the system, ju is the chemical potential for fermion number, 
and is the number of fermions. The Fermi surface is defined by a Fermi energy Ep = fj., at 
which the free energy is minimized, so adding or subtracting a single particle costs zero free 
energy. Now switch on a weak atttactive interaction. As we have just seen, it costs negligible 
free energy to add a pair of particles (or holes) close to the Fermi surface, and if they have the 
right quantum numbers then the attractive interaction between them will lower the free energy of 
the system. Many such pairs will therefore be created in the modes near the Fermi surface, and 
these pairs, being bosonic, will form a condensate. The ground state will be a superposition of 
states with all numbers of pairs, spontaneously breaking the fermion number symmetry. 

High-density low-temperature quark matter has exactly the right ingredients for the BCS 
mechanism to operate. Asymptotic freedom of QCD means that at sufficiently high density and 
low temperature, there is a Fermi surface of almost free quarks. And the interactions between 
quarks near the Fermi surface are certainly attractive in some channels, because quarks bind 
together to form baryons. We therefore expect quark matter that is sufiiciently cold and dense 
to generically exhibit color superconductivity. The densities at which the strong interaction be- 
comes perturbatively weak are extraordinarily high [4], so it remains an open question whether 
color superconducting phases persist down to the densities achieved in neutton star cores. 

The phase structure of cold quark matter is expected to be compUcated, with many competing 
phases (see Fig. 1). This is because quarks, unlike electrons, have color and flavor as well 
as spin degrees of freedom, so many different patterns of pairing are possible. Since pairs of 
quarks cannot be color singlets, the resulting condensate will break the local color symmetry 
S f/(3)coior- We therefore call Cooper pairing of quarks "color superconductivity". Note that 
the quark pairs play the same role here as the Higgs particle does in the standard model: the 
color-superconducting phase can be thought of as the Higgs phase of QCD. 

The wavefunction of a Cooper pair must be antisymmetric under exchange of the two fermions. 
The most attractive channel for two quarks is color antisymmetric (the color 3a), Dirac antisym- 
metric (the Lorentz scalar C75), and spatially symmetric (s-wave). This requires antisymmetry 
in the remaining label, flavor. We conclude that pairing between different flavors will be typically 
be the energetically favored option. As we will see, this turns out to be crucial to understanding 
the high-density phase structure of quark matter. 

2.2. Phase diagram of quark matter 

Fig. 1 (left panel) shows a schematic phase diagram for QCD that is consistent with what 
is currently known. Along the horizontal axis the temperature is zero, and the baryon density 
is zero up to the onset transition where it jumps to nuclear density; the density then rises with 
increasing //. Neutron stars are in this region of the phase diagram, although it is not known 
whether their cores are dense enough to reach the quark matter phase. Along the vertical axis 
the temperature rises, taking us through the crossover from a hadronic gas to the quark-gluon 
plasma. This is the regime explored by high-energy heavy-ion colliders. 

At the highest densities we find the color-flavor locked (CFL) color-superconducting phase, 
in which the strange quark participates symmetrically with the up and down quarks in Cooper 
pairing. The CFL phase may extend all the way down to a few times nuclear density, or there 
may, as shown in the figure, be an interval of some other phase or phases. These may include 
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Figure 1: Left panel: conjectured phase diagram of matter at ultra-high density and temperature. The CFL phase is 
a color superconductor that is superfluid (like cold nuclear matter) and has broken chiral symmetry (like the hadronic 
phase). Right panel: calculated phase diagram using an NJL model and considering only spatially uniform phases [12]. 



two-flavor color superconductivity ("2SC"), crystalline color superconductivity ("LOFF") [5, 6] 
or some form of single-flavor pairing [7, 8, 9, 10, 11]. 

2.3. Color-flavor locking ( CFL) 

At the highest densities, the favored pairing pattern is "color-flavor locking" (CFL) [13]. This 
has been confirmed by weak-coupling QCD calculations [14, 15], Dyson-Schwinger calculations 
[16], and in Nambu-Jona-Lasinio models [13, 17, 18]. The CFL pairing pattern is 

(qi<fj)iP, ^ Cy,{(K + 1)6'^^. + {K- l)5J(5f) 

[S t/(3)eolor] X S f/(3)z, X S U{2>)r XU{1)b ^ S f/(3)c+L+fi XZ2 



(1) 



3 [U{1)q\ 
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Color indices a,fi and flavor indices i,j run from 1 to 3, Dirac indices are suppressed, and C 
is the Dirac charge-conjugation matrix. The term multiplied by k corresponds to pairing in the 
(6s , 6s ), which although not energetically favored breaks no additional symmetries and so k 
is in general small but not zero [13, 14, 15, 19]. The Kronecker deltas connect color indices 
with flavor indices, so that the condensate is not invariant under color rotations, nor under flavor 
rotations, but only under simultaneous, equal and opposite, color and flavor rotations. Since color 
is only a vector symmetry, this condensate is only invariant under vector flavor-)-color rotations, 
and breaks chiral symmetry. The features of the CFL pattern of condensation are 

— The color gauge group is completely broken. All eight gluons become massive. This 
ensures that there are no infrared divergences associated with gluon propagators, so at 
asymptotically high densities this phase can be rigorously studied in perturbation theory. 

— All the quark modes are gapped. The nine quasiquarks (three colors times three flavors) 
fall into an 8 ffi 1 of the unbroken global S U{3), so there are two gap parameters. The 
singlet has a larger gap than the octet. 
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Figure 2: (Color online) Illustration of the splitting apart of the Fermi momenta of the various colors and flavors of quai'ks 
(exaggerated for easy visibility). In the unpaired phase, requirements of neutrahty and weak interaction equilibration 
cause separation of the Fermi momenta of the various flavors. The splittings increase with decreasing density, as /j 
decreases and Mj(/j) increases. At very high density the splitting is small, favoring the the CFL phase, where all colors 
and flavors pair and have a common Fermi momentum. At intermediate density we expect complicated compromises 
between pairing and Fermi-surface splitting, for example the 2SC phase, where up and down quarks of two colors pair, 
locking their Fermi momenta together. 



— A "rotated electromagnetism" survives unbroken. Its generator is Q, a linear combination 
of a color rotation and an electromagnetic phase rotation; its gauge boson is therefore a 
combination of the original photon and one of the gluons. The CFL phase is electrically 
neutral without any electrons [20], and is therefore a transparent insulator. 

— Two global symmetries are broken, the chiral symmetry and baryon number, so there are 
two gauge-invariant order parameters that distinguish the CFL phase from the QGP, and 
corresponding Goldstone bosons which are long-wavelength disturbances of the order pa- 
rameter When the light quark mass is non-zero it explicitly breaks the chiral symmetry 
and gives a mass to the chiral Goldstone octet, but the CFL phase is still a superfluid, 
distinguished by its baryon number breaking. 

— The symmetries of the 3-flavor CFL phase are the same as those one might expect for 
3-flavor hypernuclear matter [17, 21], so it is possible that there is no phase transition 
between them. 



3. Cooper pairing in the real world: 2+1 flavors 

In a real compact star we must require electromagnetic and color neutrality [22, 10, 23, 24], 
allow for equilibration under the weak interaction, and include a realistic mass for the strange 
quark. These factors tend to pull apart the Fermi momenta of the different quark species, impos- 
ing an energy cost on cross-species pairing, and hence disfavoring the CFL phase at sufficiently 
low densities. To see how this occurs, consider the left panel of Fig. 2, which shows the Fermi 
momenta of the diff'erent colors and flavors of the quark species. The strange quarks have a lower 
Fermi momentum because they are heavier, and hence are more energetically costly. To maintain 
electrical neutrality, the number of down quarks must be correspondingly increased, so the down 
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quark Fermi momentum is raised. To lowest order in the strange quark mass, the separation 
between the Fermi momenta is Spp = M^/(4/i), so the spHtting is smaller at higher densities. 
Electrons are also present in weak equilibrium, with /j.^ = M^JiAjx), so their charge density is 
parametric ally of order /z^ ~ M^/fi^ <K fiMj, meaning that they are unimportant in maintaining 
neutrality. 

In the CFL phase the situation is completely different. All the colors and flavors pair with 
each other, locking all their Fermi momenta together at a common value (Fig. 2, right panel). This 
is possible as long as the energy cost of forcing all species to have the same Fermi momentum is 
compensated by the pairing energy that is released by the formation of the Cooper pairs. Working 
to lowest order in M^, we can say that parametrically the cost is i^dp^p ~ M,, and the pairing 
energy is jj}/^^^, so we expect CFL pairing to become disfavored when Acfl ^ M^l^; actually 
the CFL phase remains favored over the unpaired phase as long as Acfl > M^JA/j [10], but 
already becomes unstable against unpairing when Acfl ^ M^/2fi [25, 26]. Sch winger-Dyson 
calculations [27] confirm that the CFL tends to be favored over other phases like the 2SC phase, 
and NJL model calculations [26, 28, 29, 30, 31] find that if the attractive interaction were strong 
enough to induce a 100 MeV CFL gap when - then the CFL phase would survive all the 
way down to the transition to nuclear matter. Otherwise, there must be a transition to some other 
quark matter phase: this is the "non-CFL" region shown schematically in Fig. 1 . 

When the stress is small, the CFL pairing can bend rather than break, developing a condensate 
of K'^ mesons, [32]. When the stress is larger, however, CFL pairing becomes disfavored. A 
comprehensive survey of possible BCS pairing patterns shows that all of them sufl'er from the 
stress of Fermi surface splitting [33], so in the intermediate-density "non-CFL" region we expect 
more exotic non-BCS pairing patterns. 

4. Compact star phenomenology 

The high density and relatively low temperature required to produce color superconducting 
quark matter may be attained in compact stars (neutron stars). This opens up the possibility of 
using astronomical observations to obtain data on the phase diagram of quark matter, although 
it must be admitted that a neutron star is not an ideal laboratory. Most of them are thousands 
of light-years from earth, and this limits the features that we can observe. Even so, there is an 
ongoing effort to develop signatures for the presence of quark matter in neutron stars (for a longer 
review see Many of these exploit the expected color superconductivity of quark matter, 
which has a a profound effect on transport properties such as mean free paths, conductivities and 
viscosities. In this section we give a brief summary, concentrating on neutron stars with quark 
matter cores, ("hybrid stars"). Pure quark matter stars ("strange stars") only exist if quark matter 
is more stable than nuclear matter even at zero pressure, and we will not discuss that possibility. 

4.1. Quark matter and the mass-radius relation 

In principle one might think that color superconductivity should affect the mass-radius re- 
lation for neutron stars with quark matter cores, ("hybrid stars") since it affects the equation of 
state (EoS) at order (A/fi)^ [34, 35]. However, other parameters such as the effective strange 
quark mass can have similar effects on the EoS, so it is hard to distinguish color-superconducting 
quark matter from unpaired quark matter using the M(R) curve. 

Actually, the M(R) curve does not clearly tell us whether there is any kind of quark matter 
in the star. Some authors have relied on the idea that quark matter is "soft", which would mean 
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that a hybrid star has a low maximum mass, and so finding a neutron star with a mass of order 
1.8 Mq or higher would rule out the presence of quark matter in its core (e.g. [36]). This is 
true for matter consisting of free quarks, but when one includes reasonable estimates of strong 
interaction corrections the quark matter EoS becomes considerably stifl^er. Hybrid stars can then 
have masses up to 2 Mq, and their M(R) curves become almost indistinguishable from those 
predicted by commonly-used models of nuclear matter [37, 38]. 

4.2. r-mode spindown 

The r-mode is a bulk flow in a rotating star that, if the bulk and shear viscosities are low 
enough, spontaneously arises and radiates away energy and angular momentum in the form of 
gravitational waves [39, 40]. Since viscosity is a sensitive function of temperature, this leads to 
"forbidden regions" in the Q.-T (spin frequency vs temperature) plane: any star that started off in 
such a region would quickly spin down via the excitation of r-modes, and exit the region. Any 
hypothesis about the interior constitution of a neutron star will lead to predictions of its viscosity, 
and hence a characteristic forbidden region in the Q.-T plane. This is illustrated in Fig. 3, which 
is taken from Ref. [41]. We see that the forbidden region for various models of nuclear matter 
(left panel) is quite different from that for models of hybrid stars (right panel). The analysis 
neglects potentially important features, such as mutual friction (phonon-vortex scattering) [42] 
and modification of the r-mode profile by the non-uniformity of the star, but illustrates how 
astrophysical observations can probe neutron star interiors. 

4.3. Quark core density discontinuity and gravitational waves 

The interface between a quark matter core and a nuclear matter mantle could be a sharp inter- 
face with a jump in energy density. (The alternative is a mixed phase with a smooth density gra- 
dient, but this only occurs if the surface tension of the interface is less than about 40 MeV/fm^ - 
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0.2 X (200 MeV)-', a fairly small value compared to the relevant scales Aqcd ~ 200 MeV, 
H ~ 400 MeV [43].) A sharp interface might modify the signature of gravitational waves emitted 
during mergers and detected via observatories like LIGO, since those encode information about 
the ratio MIR [44], and the star would in effect have two radii, one for the quark core and one for 
the whole star. 

4.4. Crystalline pairing, gravitational waves, and pulsar glitches 

One candidate for the intermediate "non-CFL" quark matter phase of Fig. 1 is the "LOFF" 
crystalline phase [5]. Current indications are that the crystal has a much higher shear modulus 
(v ~ 0.5-20 MeV/fm^) than nuclear matter (v ~ 10 ^^-lO'^ MeV/fm^) [45]. One resultant 
signature is that the quark matter is rigid enough to sustain a large quadrupole moment, leading to 
detectable emission of gravitational waves. The LIGO non-detection of such gravity waves from 
nearby neutron stars already shows that they do not have quark matter cores that are deformed to 
the maximum extent allowed by the estimated shear modulus [46, 47]. 

Two other relevant phenomena are glitches, in which pulsars speed up their rotation occa- 
sionally, and precession. However, it is hard to come up with a mechanism that allows for both 
these phenomena in the same star [48]. The standard glitch mechanism involves piiming of su- 
perfluid vortices in the crust, which would suppress precession in all stars, since they all have 
crusts. Quark matter offers a way out — glitches could arise from pinning in a crystalline quark 
matter core. Then there would be two populations: heavy stars with a crystalline core which 
could glitch but not precess; and lighter stars with no core which could precess but not glitch. To 
test this we need better calculations of the properties of the crystalline phase and more detailed 
observations of glitch rates and precession frequencies. 

4.5. Cooling by neutrino emission 

The cooUng rate is determined by the heat capacity and emissivity, both of which are sensitive 
to the spectrum of low-energy excitations, and hence to color superconductivity. CFL quark 
matter, where all modes are gapped, has a much smaller neutrino emissivity and heat capacity 
than nuclear matter, and hence the cooling of a compact star is Ukely to be dominated by the 
nuclear mantle rather than the CFL core [49, 50]. Other phases such as 2SC or LOFF give 
large gaps to only some of the quarks. Their cooling would proceed quickly, then slow down 
suddenly when the temperature fell below the smallest of the small weak-channel gaps. This 
behavior should be observable [51, 52, 53]. There is already evidence that, although the cooling 
of many neutron stars is broadly consistent with the standard cooling curves, some fraction of 
neutron stars cool much more quickly [54]. One may speculate that Ughter neutron stars cool 
following the standard cooling curve and are composed of nuclear matter throughout, whereas 
heavier neutron stars cool faster because they contain some form of dense matter that can radiate 
neutrinos via the direct Urea process [55]. This could be quark matter in one of the non-CFL 
color-superconducting phases, but there are other, baryonic, possibilities. If this speculation is 
correct, then if neutron stars contain CFL cores they must be "inner cores", within an outer core 
made of whatever is responsible for the rapid neutrino emission. 
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